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rn '. Abstract 

The Super-Kamiokande best global fit, which includes data from SNO, Gallium and Chlorine 

> ; 

■ experiments, results in a hep neutrino contribution to the signals that, even after oscillation, is 

o 

(N 



greater than the SSM prediction. The solar hep neutrino flux that would yield this contribution 



\ is four times larger than the one predicted by the SSM. Recent detailed calculations exclude that 

o ■ 

the astrophysical factor 5/iep(0) could be wrong by such a large factor. Given the reliability of the 



temperature and densities profiles inside the Sun, this experimental result indicates that plasma 
1) , effects are important for this reaction. We show that a slight enhancement of the high-energy 

tail, enhancement that is of the order of the deviations from the Maxwell-Boltzmann distribution 



^ ■ expected in the solar core plasma, produces an increment of the hep rate of the magnitude required. 

We verified that the other neutrino fluxes remain compatible with experimental signals and SSM 
predictions. Better measurements of the high-energy tail of the neutrino spectrum would improve 
our understanding of reaction rates in the solar plasma. 
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I. INTRODUCTION 



Results of the SK and SNO experiments JI], provide model-independent evidence of 
electon-neutrino oscillation into other flavors. A recent global fit of the Super-Kamiokande 
(SK) Collaboration 0] that includes data of SNO, Gallex/GNO and SAGE shows that 
the Large Mixing Angle (LMA) solution is preferred at the 98,9% confidence level. The 
corresponding mass-square difference Am^ is within the range 3 ■ lO-^eV^ < Am2 < 19 • 
10~^eV^, and the mixing angle 9 within the range 0.25 < tan^ 9 < 0.65. 

The flux resulting from the SK best fit, (5.33 ± 0.36) ■ 10^ cm^^s"^, is in substantial 
agreement with the SSM value (5.05lo8i) ' 10^ cm~^s~^. Also the flux measured by 
SNO 1^ assuming the standard ^B energy spectrum, (5.091q;43) ' cm~^s~^, or the one 
obtained using a distorted spectrum, (6.42 ± 1.57) ■ 10^ cm~^s~^, are in agreement with the 
SSM value and with the SK fit. On the contrary, the hep flux resulting from the SK best 
fit, 

= 36 ■ 10 W , 

which corresponds to the LMA solution, is about four times higher than the one predicted 
in the SSM, 9.3 • 10^ [i, if one uses ShepiO) = 10.1 ■ lO^^o keV barn §, or 

= 7.9 • 10^ cm- , 

if one uses the more recent Shepi^) = (8.6 ± 1.3) ■ 10^^° keV barn |^. The second best 
fit (LOW) gives also a hep fiux considerable larger than the one predicted in the SSM, 
23 ■ 10^ cm^^s^^. Only the Quasi- VAC and SMA solutions, which appear disfavored from 
present data, give a hep fiux compatible with the SSM 0. 

Given the actual experimental evidence, we should justify a hep fiux four times larger 
than in the SSM. Bahcall and Krastev discussing the possibility that the hep fiux 
be considerable larger than in the SSM, conclude that they cannot find a first-principle 
physical argument demonstrating that the astrophysical factor Shep{0) cannot exceed 20 
or 40 ■ 10"^° keV barn, values which could explain the experimental value of the hep fiux. 
However, the latest detailed calculations H, ^ exclude an astrophysical factor of this order 
of magnitude; the latest result gives 

ShepiO) = (8.6 ± 1.3) ■ 10"^° keV barn . 
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Excluding that the cross section be wrong by a factor of four, the average rate can be 
four times larger only if the product of the thermal average of the cross section {av) times 
the p and ^He densities is four times larger. But both the temperature and the p density are 
accurately measuread by helioseismology in the region where hep neutrino are pro- 

duced. In addition increasing the temperature or the densities would affect much more the 
^Be-^ ^B- and CNO-neutrino fluxes IjTTj. There remains the possibility of some improbable 
mechanism that increases the ^He density, which is not measured by helioseismology [1T2| , 
by a factor of four in the /lep- neutrino production region {t/Rq > 0.12) but not in the ''-Be-, 
and ^B-neutrino production regions {t/Rq < 0.12). Note that mixing would produce the 
opposite effect [|l^ . 

We propose that the large increase of hep flux which flts the present experiemental data is 
a consequence of the enhancement the high-energy tail of the ^He-p momentum distribution. 
The required hep flux of 36 ■ 10^ cm~^s~^ is a signal that the thermal ion distribution of 
the solar plasma deviates slightly from the standard one in the region where the reaction 
^He + p —>-^ He + e"*" + z/g is active. In fact small deviations in the high-energy tail give 
dramatic effects on the fusion rates without affecting the bulk properties of the medium, 



which are the ones measured by helioseismology |T3| , |14|. 

Let us remind a few of the reasons that lead to deviations from the Maxwellian distri- 
bution [^, In plasmas with parameters such as those that can be found in part of 
the solar core, in the solar atmosphere or in the interior of giant planets |T^, the Debye 
screening is approximately valid, but the essential many-body character of the interaction 
must be taken into account (nonlocality); the time necessary to build up again screening 
after hard collisions is comparable to the inverse of plasma frequency (memory effect or 
time nonlocality); since many collisions are necessary before particles loose memory of their 
initial state, the scattering process cannot be considered Markovian. The conditions are 
such that there is no clear separation between collective and individual degrees of freedom 
and the description of the system in the relevant energy range requires the use of additional 
scales, dynamically generated by the interaction, in addition to the temperature: the pure 
exponential exp{—E/kT), which is determined by the single scale kT, is insufficient. The 
Maxwellian regime ends below a given scale. It has also be shown ||T6[ that a flnite width of 
the quasi-particle is effectively equivalent to a distortion of the momentum distribution. 

Small deviations from the Maxwellian distribution, as the ones that are relevant for 
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us, can be parameterized by the deformation parameter S introduced many years ago by 
Clayton et al. ||T^. This description could be related to the general framework for non- 



Maxwellian distributions that has been proposed by Tsallis and applied in many different 
contexts [113, 113: the corresponding distribution depends on parameter q that can be related 



to 5 in the appropriate limit. The parameter 5 can be related to the plasma parameter F 
and the ion-ion correlation special parameter a |]IB|, |T^|, and to the finite width of the 
quasi-particle [|16 |. 



In the next Section we derive the hep flux. In the third Section we present the expression 
that relates the deformation parameter 5 or g to the plasma parameter and the ion-ion cor- 
relation special parameter a; in addition we calculate the effect on the other neutrino fluxes, 
considering the LMA oscillation and the admissible deviations from Maxwellian distribution 
of the different reacting ions. In the last Section we summarize and comment our results. 

II. THE hep FLUX 

The reaction ^He + p —^^ He + e^ + Vf. produces a neutrino with an endpoint energy of 
18.81 MeV, the highest energy expected for solar neutrinos and the only one above 15 MeV. 
The rate of the hep reaction is very slow and does not affect the solar structure. 

The SSM hep flux, which uses standard distribution in the plasma, is [§: 

- 21 ■ (1 ± 0.03) . ( ,3.J-4°i,^J • 10» cn-.- . 



Latest values of ShepiS^) have ranged from (2.3 ± 0.9) ■ 10"^" keV barn to (10.1 ± 0.9) • 
10^^° 0. Let us mention that Alberico et al. pT| have proposed a method to determine 



'S'hep(O) from a (challanging) laboratory experiment on electron scattering. The latest and 
more reliable determination 

Shepi'd) = (8.6 ± 1.3) • 10-^\eV barn 

gives 

= (7.9 ± 1.4) ■ 10^ cm-^s-i , 

where the quoted error includes only the contribution due to Shepi^). 

Super-Kamiokande best global fit (LMA) yields a hep flux, 36 ■ 10^ cm~^ s~^, about four 
times $}Jep ' while their second best fit (only 1% of probability) yields 23 ■ lO'^ cm~^ s~\ 



about three times $^ep |T[; therefore 



C'p = (3^4)xCp • 

Fusion rates are proportional to the local densities of incoming particles and to the 
average cross section {va{v)) where v is the relative velocity. Therefore, rates depend on 
the statistical distribution, i.e., on the local temperature if the distribution is Maxwellian, 
and also on at least an additional deformation parameter if the distribution deviates from 
the Maxwellian. For example, small deviations can be parameterized by the form proposed 
by Clayton 0: 

^-E/kT _^ ^-E/kT-S{E/kTf (^-^^ 

for 5 > 0; this distribution (Druyvenstein distribution) introduces a scale kT/5 in addition 
to kT\ the second term in the exponential becomes important when E > kT/5. 

To first approximation, we disregard the changes of the other rates and write the change 
of the average hep cross section, and therefore of the rate, as [|l^, 



{va{v)),^^^ = {va{v)),e-^^'^'-'^ (2) 

where 71,3 = {Eo/kT)l^p with Eq = {EaikTf /Af'^; Eq = 2fxc'^{ZiZ2a7rf is the Gamow 
energy and {va{v))Q the Maxwellian rate. The indices (1,3) indicate the reaction + 
^He. For this reaction 71 3 = (9.07 keV/1.01 keV)^ ~ 81. For small enhancement of the tail 
of the distribution the deformation parameter S is negative and related to Tsallis entropic 
parameter q, 6 = {1 — q)/2 with g > 1. Note that even if the distribution of Eq. (|lD is 
normalized only for 5 > 0, the resulting change of the rate in Eq. (0), which is derived as an 
asymptotic expansion in S, can be analytically continued to 5 < 0. 
The required enhancement 

g-7i,35i,3 = 3 ^ 4 implies -0.017<(5i,3 < -0.014 . 
Let us compare this range of values of S with the relation derived and discussed in 



Refs. 13, 14 



\5ij\ = UafjTiTj , 

which gives 5ij in terms of Fj = Zi'^j Zje^{njY^^ /kT, the plasma (or ion-ion) coupling 
constants (Fj ^ 0.14Zj in the Sun core), where nj is the number density of ion j, and 
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of Ojj, the ion-sphere- model parameters introduced by Ichimaru ||22|, parameters that 
are related to the ion-ion correlation fuctions (0.4 < a < 0.9). An alternative formula 
by Ichimaru which uses F^^ = (^ZiZjC^ '^i''^) /^-^)) instead of TiFj, yields 

similar ion-sphere-model paramenters. 



III. THE OTHER FLUXES 

In this Section we want to give a quantitative estimate of the range of 5's compatible 
with the experimental determination of the main neutrino fluxes, and compare it with the 
range of values suggested for the hep neutrino flux. 

Previous calculations that assumed standard neutrinos have shown that 5 of the order of a 
few times 10^^ were compatible with the neutrino experiments even if these deviations from 
the Maxwellian distribution were not sufficient to solve the solar netrino problem [113], 11^ [16 



Since we now have experimental evidence for oscillations and the best fit to the oscillation 
parameters (LMA) yields fiuxes close to the one predicted in the SSM, we expect that the 
values of 5 compatible with present data be even closer to zero. Therefore, the effect on the 
neutrino fiuxes of small deviations from the Maxwellian distribution can be reliably obtained 
by using power laws that include the reajustment of the solar model 



We repeat the analysis of Ref. [|T4| taking into account the oscillation probability. We add 



the survival probability factors Pi of i-th neutrino fiux to the four equations of Castellani et 



al. [^]. The Gallium ||2^, Chlorine ||26| and Luminosity equations yield: 

70.8 ± 4.5 = 73.185 • RppPpp + 34.196 ■ RsePBe + 9.0165 ■ RcnqPcno + 2.43 ■ $|^^3a) 
2.56 ± 0.23 = 0.0226 ■ RppPpp + 1.1448 ■ RBePse + 0.4239 ■ RcnqPcno + l-H ■ $|^^3b) 
63.85 = 0.980 ■ 59.85 ■i?pp + 0.939- 4.77 

+0.937 ■ 1.034 ■ RcNO + 0.498 ■ 10"^ ■ 5.05 ■ Rb , (3c) 

where 5.05 is the produced fiux; the fiux in units of 10^ cm~^ s~^ 

= (5.05lJ:°J) ■ RbPb = 1.76 ± 0.05 (3d) 

is measured in the experiment SNO. The SSM fiuxes $ are the ones calculated by Bahcall 
et al. 0] . We use the survival probability Pj extracted from the experiements by Berezinsky 



and Lissia |28|: Ppp = 0.58 at 0.265 MeV, Pbc = Pcno = 0.55 at 0.814 MeV, and Pb = 0.32 



at 6.71 MeV. The quantity defined in Ref. |T3[, represents the rate enhancement or 
depletion due to the nonmaxwelhan tail. Eliminating Rpp and using the relation between R 
and 6 [0, |T^ , the other i?j (with conservative errors) are: 

Rj^^ = e235pp-139.755He = 0.87 ± 0.09 

RcNo = e62.15pp+6.55ffe-4075p = 0.87 ± 0.09 (4) 
= e62.15pp+149.55He-1905p ^ 1.10 ±0.20 

where for the purpose of this estimation we have introduced several 5's: 6pp for protons 
in the pp production region {t/Rq > 1), which is also the region where the hep neutrinos 
are produced, 5p for protons in the inner core {t/Rq < 1), where all other neutrinos are 
produced, and Shs for the distribution of the He nuclei ('^H and ^H). The result is that 
6pp ~ -10-2, while \5p\ ~ \5He\ < 2 ■ IQ-^. 

We note that neutrinos from the pp and hep reactions are produced in the same region 
(r/i?0 > 1), and they have similar values (limit) for 6. The difference is that values of 
6 ~ 2% have large effect on the hep reaction |]16|, while they have a much smaller effect on 
the pp reaction both because the Coulomb barrier is lower and because the adjustment of 
the solar structure to a change of the pp rate, which is directly linked to the luminosity, is 
such that it minimizes the change of the rate, that ends up being between 0.5% and 3.0% 
larger that in the SSM (at the solar surface). 

The other neutrino fluxes are all produced in the inner core {t/Rq < 1), where the 
plasma conditions are different, and they show much smaller deviation from Maxwellian 
distribution. We note that the value of a deduced from the 6 of the proton component, 
i.e, a value between 0.54 and 0.72, is in agreement with the range of values allowed for a 
weakly coupled and weakly non ideal plasma (F < 1). In addition the values of the different 
parameters 6 are within the constraints imposed by helioseismology, as it was checked in 



Ref. [£9| without considering oscillations: oscillations makes the agreement even better. In 
detail protons superdiffusion in the inner core would yield a ^Be flux reduced by 4-22% with 
respect to SSM calculations; the same range of reduction is obtained for the CNO flux (these 
fluxes are at the solar surface, before oscillation). 

Such values of 6 could be measured or excluded by more precise measurements and solar 
model calculations. 
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IV. DISCUSSION AND CONCLUSION 



Our work is based on the following points: 

• the global best fit (LMA) to the neutrino experiments requires a /lep-neutrino produced 
flux of 36-10^ cm~^ s^^ |jl||, about four times the one predicted by the SSM also the 
second best fit (LOW) requires a flux three times the one in the SSM; if one consider 
the contribution to the signals on earth from hep neutrinos, this contribution, even 
after oscillation is larger than the SSM prediction; 

• the latest precise evaluation [0 of the astrophysical factor, Shep{0) = (8.6 ± 1.3) • 
10~^° keV barn, which is slighter smaller than the values used in the SSM 0, rules 
out that the cross section could be a factor of four larger; 

• temperature and p density are measured very precisely by helioseismology in the region 
of interest and could not vary significantly without dramatically affecting the ^B- and 
^Be-neutrino fluxes; 

• it is not known a physical mechanism that could increase the density of '^He in the 
/lep-neutrino production region without increasing the same density in the ^B- and 
''Be-neutrino production region (any mixing would produce the opposite effect) [p!2| ; 

• there exist several mechanisms that deform the high-energy tail of the ion velocity 
distribution in plasmas ||13|, 0, |T6|, ^; we do not exclude that other mechanisms 



could be found that also produce deformation of the distribution. 

We have proposed that the experimental results are signals of plasma effects on the fusion 
rates. 

We have verified that the slight enhancement of the high-energy tail of the distribution 
naturally produces the needed large /lep-neutrino flux; the resulting deformation parameter 
is in agreement with the theory of weakly coupled and weakly non-ideal plasmas. This 
kind of modifications of the distribution do not affect the bulk properties and, in particular, 
the measured helioseismiologic quantities. The other neutrino fluxes are not significantly 
modified. 

Future more precise experimental determinations of the oscillation parameters, Am^ and 
mixing angle 6, and of the high-energy neutrino spectrum (and in particular a direct mea- 
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surement of the hep contribution) and improved theoretical calculations, or experimental 
determination |21|, of the astrophysical factor Shep{0) could corroborate (or bound) the 
importance of plasma effects for some of the fusion rates in the Sun. 
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